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ABSTRACT 

We derive simple analytic expressions for the flux and spectrum of ultra-high energy cosmic -rays (UHECRs) 
predicted in models where the CRs are protons produced by extra-Galactic sources. For a power-law scaling 
of the CR production rate with redshift and energy, dri/de oc e~"(l our results are accurate at high en- 
ergy, e > 10'^^ eV, to better than 15%, providing a simple and straightforward method for inferring dri/deie) 
from the observed flux at e. We show that current measurements of the UHECR spectrum, including the latest 
Auger data, imply e^dn/deiz = 0) = (0.45 ±0.15)(a- 1) x lO'*'* erg Mpc"^ yr"' at e > 10''^'^eV with a roughly 
confined to 2 < a < 2.7. The uncertainty is dominated by the systematic and statistic errors in the experimental 
determination of individual CR event energy, (Ae/e)sys ^ (Ae/£)stat ~ 20%. At lower energy, dri/de is uncer- 
tain due to the unknown Galactic contribution. Simple models in which a ~ 2 and the transition from Galactic 
to extra-Galactic sources takes place at the "ankle", e ^ 10'''eV, are consistent with the data. Models in which 
the transition occurs at lower energies require a high degree of fine tuning and a steep spectrum, a ~ 2.7, 
which is disfavored by the data. We point out that in the absence of accurate composition measurements, the 
(all particle) energy spectrum alone cannot be used to infer the detailed spectral shapes of the Galactic and 
extra-Galactic contributions. 
Subject headings: cosmic rays 



1. INTRODUCTION 

The origin of the observed Cosmic Rays (CRs) at differ- 
ent ene rgies is st ill unknown (se e Blandford & Eichler 1987; 
lAxford 1994; N agano & Watsonll2000l for reviews). The cos- 
mic ray spectrum changes its qualitative behavior as a func- 
tion of particle energy; it steepens around 5 x 10'^ eV (the 
"knee") and flattens around 5 x 10'** eV (the "ankle"). Be- 
low ^ 10'^ eV, the cosmic rays are thought to originate from 
Galactic supernovae. The composition is dominated by pro- 
tons at the lowest energies, and the fraction of heavy nuclei 
increases with en ergy. The proton fraction at ~ 10^^ ey is re - 
duced to - 15% (iBrunett et. allll990l:lBemlohr et. allll998l) . 
At yet higher energies, there is evidence that the fraction 
of light nuclei increases, and that the cosmic -ray flux above 
5 X 10'^ eV is again dominated by protons (Gaisser et. al 
1998; Bird et al. 1994; note that preliminary results by the 
Auger collaboration suggest that at the highest energies the 
composition becomes heavier again, e.g. Bluemer; for the 
Pierre Auger Collaboration 2008). The composition change 
and the flattening of the spectrum around 10'^ eV (see Fig. 
|9|l suggest that the flux above and below this energy is domi- 
nated by different sources. At energies of £i9 = e/lO'^eV 1 
the Larmor radius of CRs in the Galactic magnetic field is 



ergy regime lO'^eV < e < 3 x 10'**eV (lLeutholdll2008l and 
references therein), in conflict with earlier results. 

CRs with energies exceeding e>5xlO'^eVare termed 
Ultra High Energy Cosmic Rays (UHECRs). The (probably 
XG) sources of these particles are unknown and have been 
the is sue of much debate (e.g. iWaxmanI 12004) ; iBerezinskvl 
l2008l and references therein). Measurements of the fluxes 
of particles in this energy range ha ve been conducte d 
by several groups including AGASA (Takeda et al.' '1998"), 
Fly's Eye (Bird et al. 1994), HiRes (Abbasi et al. 2008), 
Yakutsk (j^ fanasiev et al. 19 93i) and most recen tly Auger 
jBluemer; for the Pierre Auger Collaboratioiill2008h . 

The problem of inferring the generation rate and generation 
spectrum of UHECRs from the observed flux and spectrum 
has been addressed by several authors (e.g. [V v^axman 19951 
'Bahcall& Waxman"2003l: iBerezinskv et al.ll20b6t IBerezinskvl 
i2008; Stanev 2008, and references therein). As long as we 
restrict to the highest UHECRs the di fferent models ro ughly 
agree with each other For example, IWaxmanI ( Il995h finds 
a generation spectrum dti/de = Q{e) oc e~" with a = -2.3 ± 
0.5 and energy production rate of CRs in the energy range 
lO'^^eV < e < lO^'eV of 4.5 ± 1 5 x 10^** erg Mp c "^ yr" ' 



• 3B_5 gEigZ ' kpc, 



(1) 



where B = 3B-5.5/iG is the value of the Galactic magnetic 
field and Z is the charge of the observed Galactic cosmic 
rays. Since the Galactic magnetic field can not confine pro- 
tons above 10'^ eV, it is believed that the nearly isotropic 
cosmic ray flux at e > 5 x 10''^ eV originates from extra- 
Galactic (XG) sources. The small, but statistically significant, 
enhance ment of the flux at e < 3 x lO'^ eV near the Galac- 
tic plane (iBird et al.i ri999; Hay shida et al. l fT999). suggests a 
Galactic origin at these lower energies. Note, however, that 
the Auger experiment did not detect any anisotropy in the en- 
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Israel; 



while in the model presented by IBerezinskv et alj (l2006l) . 
the spectral index is a = 2.1 and the energy generation rate 
in the energies 10''^eV < e < 10^'eV is approximately 6 x 
10'*'*ergMpc-^ yr-i. The main difference between differ- 
ent models is the location of the transition between Galac- 
tic dominated and XG dominated CRs. While some models 
assume that the transition is at energies of e ~ lO'^eV (e.g. 
Wax man 1995), other models alow for a transition at con- 
sidera bly lower energies, e < lO'^eV (e.g. IBerezinskv et alJ 
l2006h . 

In this paper we revisit this problem, presenting a novel an- 
alytic tool for the analysis of the UHECR spectrum and in- 
cluding the most recent data from the Auger experiment. Our 
estimates for the energy generation rate, are in agreement with 
previous estimates (see §[3]l. 
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This paper is organized as follows. We first describe the 
framework of our analysis in § |2] We focus on models where 
the UHECRs are dominated by protons originating from XG 
sources with power-law generation spectrum and an energy 
production rate that depends on redshift. We derive simple 
analytic approximations for the calculation of the effects of 
the energy losses due to interaction with the CMB, and show 
that the approximate analytic results reproduce with good ac- 
curacy the results of direct numerical calculations for particle 
energies of e > 10'**'^eV. Next, we estimate in §[3] the en- 
ergy production rate of cosmic rays with energies e > lO'^eV, 
and show that it is roughly model independent. In section 
§ |4] we show that the measured spectrum of CRs at energies 
e > lO'^eV is consistent with simple models with a transition 
between Galactic and XG dominated regions at the 'ankle'. 
We review and illustrate the argument that models in which 
the transition between Galactic and XG sources occurs at en- 
ergies considerably below the ankle require a high degree of 
fine tuning. We summarize the results and the main conclu- 
sions in § ID 

2. ANALYTIC CR PROPAGATION MODEL 

In this section we derive approximate analytic expressions 
for the expected UHECR energy dependent flux assuming it is 
dominated by XG proton sources. First we derive in § l2.1l ap- 
proximate expressions for the flux in terms of the energy loss 
time T by assuming that at high energies r is much shorter 
than the Hubble time H^^ . The approximate expressions are 
derived by calculating the effect of the losses to first order 
in the small parameter Hqt. We then obtain in § 12.21 a phe- 
nomenological analytic approximation for the loss time of 
protons due to pair production and pion production and use 
it to derive an analytic expression for the expected flux. We 
show that these expressions are in good agreement with de- 
tailed numerical calculations. 

2.1. Analytic estimate of the effects of CR energy losses on 
the observed spectrum 

As they propagate, high-energy protons lose energy as a re- 
sult of the cosmological redshift and as a result of production 
of pions and e + e- pairs in interactions with CMB photons. 
Here we approximate the energy loss caused by scattering, 
which is a random process, as a continuous energy loss (CEL). 
In this approximation, the energy loss rate of a proton of given 
energy is taken to be the mean loss rate of an ensemble of 
protons of the same energy. The CEL approximation is excel- 
lent for pair production, for which the mean free path is small 
and the relative energy loss in a single scattering, of order 
nie/nip ^ 10"^, is also small. For pion production, the average 
relative energy loss in a single collision is 0.13 at the thresh- 
old and rises to 0.5 at higher energy. However, fluctuations in 
proton energy resulting from this process are significant only 
for propagation dista nces smaller than 100 Mpc and proton 
energies > 10-" eV (lAharonian & CroninlfT994<) . The CEL 
approximation therefore gives accurate results for the flux be- 
low 10^" eV. At higher energies, the flux obtained by this ap- 
proximation drops faster than that obtained when fluctuations 
in proton energy are taken into account. It has been s hown 
jBerezinskv et al.ll 19751; iBerezinskii & Grigor'e"valll988 ) that 
for a flat generation spectrum, a < 2.6, the flux obtained us- 
ing the CEL approximation is accurate to better than 10% up 
to 3 X 10^° eV, the highest energy at which events have been 
reported. 



For numerical estimates we assume the following cosmo- 
logical parameters: Ho = 72 km sec"' Mpc~', VIm = 0.28 and 
J7i = 0.72. We assume that the generation spectrum is a power 
law and that the production rate depends on z as 

dri/ds = Q(e,z) = (£/£())""(! +z)"'Q{eo) (2) 

for z < Zmax and e < e,„f„, where Zmax and represent the 
maximal z and particle energy respectively, up to which CRs 
are generated. In all numerical calculations we adopt z,mx = 4 
and Emax = lO^^eV. For energies e > lO'^'^eV the results are 
insensitive to these choices. In all numerical calculations we 
use the loss time calculated bv Berezinskv et al. (2006). 

Consider a CR that propagated through the IGM while suf- 
fering energy losses due to interactions with the CMB and due 
to the expansion of the universe and reached us at f = with an 
energy eq. Denote the energy that it had at an earlier time t by 
e'(eo,t). The time dependent energy e' satisfies the equation: 



1 c/z^ 



(3) 



where TcMB(£',t) is the energy loss time due to interactions 
with the CMB and the second term on the rhs is the energy 
loss rate due to the Hubble expansion. The spectrum of parti- 
cles observed at the current time is given by: 

de' 

dtQ(e'{e(),t),t)—\,deo (4) 



^(eo)deQ ■■ 
de 



deo 

de' de' 

-Q(e'(eo,t),t)—\,deQ, 



e'T-\e',t)' 



deo 



(5) 



where Q is the CR generation rate per comoving volume per 
unit energy. 

At high energies the loss time due to the interactions with 
the CMB is much shorter than the Hubble time: tcmb ^ ' ■ 
In what follows, we use the assumption that HoTqmb ^ 1 to 
find an approximate relation between the observed spectrum 
and the CR production spectrum. 

The CMB spectrum scales with redshift as «cMB(e,z) = (1 + 
z)'^«cmb(£/(1 and thus all continuous energy losses due 

to interactions with the CMB have the following scaling, 

TcMBie.z) = (1 +z)"VcMs((l +z)e). (6) 

It is useful to consider the evolution of the scaled energy e = 
(l+z)e. We find 



1 de 



2 dz 



edt -^^"""^''^^^l+zdt 



' {l+zrTclisim+2HQTcMB(s)), 

(7) 

where we neglected second order terms in HqTcmb{£)- 

Next, we replace the time variable with a time parameter t 
that satisfies dt = ( \ +z)^dt and obtain 



e at 



(8) 



Since the time derivative of e does not depend on time, we can 
write 

dn , f°° dedt 

—(eo)deQ = dto I de-^^Q(e,t), (9) 



where 



dto = deo - 



TCMB(eo) 



l+2HoTcMBi£o) 

Since f is a function of f only, we have 

dedt _de dt _ ^j_^^-j-4 
dedi de dt 



(10) 



(11) 
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and thus 
dn 

3-(£o) ' 
ae 



1+2//oTcmb(£o) 



(l+z)-'Q( — ,z)deo. (12) 
1+z 



For a generation spectrum Q(e,z) = (e/eo) "(1 +z)'"C(eo) 
we find 



dn tcmb(£o) , 



(1+z) 



—4+m+a 



(13) 

It is useful to express the relation between the generation 
spectrum and the observed flux in terms of an effective CR 
generation time Jeff defined by 



dn 

3-(eo) = G(eo)feff(eo)- 
de 



Equation ( fT3] ) can then be written as 

tcmb(£o) r 

1+2Hotcmb(£o) Jet, 



feff(£o) ' 



(l+z) 



-4+m+a. 



(14) 



, (15) 



where we assumed that replacing £n,ax with e,„„\- in the inte- 
gration upper limit results in a small error which can be con- 
sidered second order For the special case a + m = 4 we find 
that 



feff(eo) ' 



1 



1- 



^max 

£() 



^a-lj 1+2HoTcmb(£o) 

"(16) 

For most practical purposes, assuming that £„iax is large 
enough, we can ignore the last term on the rhs of Eq. ( fT6] l 
to obtain 



feff(eo) ' 



1 



tcmb(£o) 



a- I J 1+2HoTcmb{£o) 



(17) 



We see that in case a + m = 4 the effective generation time 
depends to first order only on the value of tcmb at eo- In fig- 
ure[T]we compare the results of numerical calculations of the 
propagation of CRs to what we expect based on equation ( fT6] l 
for different values of a and m that satisfy a+m = 4. As can be 
seen in the figure, for a > 2, the differences at large energies 
e > 5 X 10'^ are smafler than ~ 15%. 

In case a + m we need to evaluate z(e, £o)- To first order 
in HqTcmb we have: 



de 

dz ~ Hodt w HqTcmb{£} — 

e 



so 



HoTcmb(£)- 



de 



(18) 



(19) 



Using equations ( fT9] l and ( fT2] i. we find that to first order in 
Hqtq (after some straight-forward algebraic manipulation) 



feff(eo) ~ [l+//orcMB(e{))(2 + (4-a-m)/r(eo,a))] 
a- 1 ■' 



where 



/r(«) = 



' , 7"CMB,o(^£o) - 

dx — —X 

tcmb{£q) 



(20) 
(21) 



Assuming that tcmb{£) is a decreasing function of energy, we 
have 

0</(a)<^. (22) 
a- 1 



1.15 




1.05 



0.95 



E(eV) 



Fig. 1 . — Comparison of equation jl6t with direct numerical integration of 
the propagation equations for a+m — 4. We use the numerical estimates for 
the energy loss time, tcmb, given bv lBerezinskv et al.l (2063) . 

The values of fr for different values of a and different ener- 
gies are presented in figure |2] (full lines). In this figure, the 
values of the upper limit 1/(q: - 1) are shown as dashed lines. 
In figure [3] we compare the results of numerical calculations 
of the propagation of CRs to the results of equation (|20| ). As 
can be seen in the figure, at large energies, e > 5 x 10'^, Eq. 
( l20l i is accurate to better than 15%. 




E(eV) 

Fig. 2. — Numerical values of /T(a,e), given by Eq. E), for energy losses 
of protons interacting with the CMB (full lines). Dashed lines are the upper 
limits l/(a-l). 



2.2. Energy losses of proton CRs due to interactions with 
the CMB 

The loss time of protons due to interactions with the CMB 
that produce pairs and pions has been calculated in several 
publications. These processes are characterized by the exis- 
tence of a threshold energy of the photons in the rest frame 
of the protons, that is required in order to create the products 
of the interaction. The typical energies required from a pro- 
ton in order to produce an electron positron pair and a pion 
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Fig. 3. — Comparison of equation i20\ with direct numerical calculations 
for a+m ^4. 



by interacting with a photon with typical energy 3Tcmb are 
roughly 



Sth. 



• ImpnieC^/OTcMB) 1 .4 x lO'^eV (23) 



and 



£f;,,7r = mpin^.c'* /OTcmb) ^2x lO^^eV 



(24) 



respectively. 

For protons with energies that are much smaller, photons 
with energies that are much larger are required in order to 
create the products. The numbers of contributing photons de- 
creases exponentially as the proton energy is decre ased. Mo- 
tivate d by this fact, we use the following anzatz (IWaxman I 
Il995h to parameterize the energy loss time: 



T ' (e) = Tolp exp(-ec,e.p/e) + '^0,1 exp(-e,.,^/£), 



(25) 



where the values of TQ gp/.^ and Sc^ep/TT are chosen to fit nu- 
merical calculations of the loss time. By setting e^^ep = 
2.7 X 10i**eV, To^ep = 3.4 X 10'' yr; e,,^ = 3.2 x lO^OeV and 
Tf),.^ = 2.2 X 10^ yr, equ ation (l25l l reproduces the results of 
iBerezinskv et alj (l2006l) to better than 10% accuracy in the 
energy range lO'^'^ < e < 10^**^ (t his a pproximation for the 
pion production agrees with that of IWaxman I d 19951) to better 
than 10%). 

The values for Ec.ep and ec,7r are of the same order of mag- 
nitude as the threshold energies in jTM and ( |23T l. To get a 
feeling for where the typical values for r^p and t-^ come from, 
we can make crude estimates as follows: For both processes 
we express the typical loss time in terms of a typical cross 
section and relative energy loss 



T ^ (ncMBO-CT] ^) 



(26) 



where for electron positron production we have roughly aep ' 

a^ar and rjgp ^ nie/mp so 



Tep ^ (ncMBaeO-Tcmp/me) 



10^ yr, 



(27) 



while for pion production we have roughly a^^ ^ iji/ m^cf- 
and Ty^r ~ rriT, /nip so 



Ttt ~ (ncMB(Ji/mAcfmp/m^) ' ~ 7 x 10^ yr. 



(28) 



We see that the crude estimates in equations ( l28T l and (l27T l. are 
of the same order of magnitude as the fitted values. 
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We can get a better fit by using more parameters. Using the 
following anzatz: 

r"'(e) = (TQ,,pe\p(ec.ep/e) + dTep)~^ + (TQ,.^s\p(ec,^/e) + dTyry 

(29) 

and setting the values Ec.ep = 9.1 x 10'**eV, tq^^p = 0.5 x 10^ yr, 
dT,p = 3 X 10'^ yr; e,,^ = 3.5 x lO^^eV, tq'^ = 1.4 x 10^ yr, 
an d dTjr = 2.4 x 10^ yr, eq uation (|29l l reproduces the results 
of IBerezinskv et al\ (l2006i) to better than 3% accuracy in the 
energy range lO'^-^ < £ < 10^° ''. 

For most purposes it is sufficient to consider only the simple 
cases where a + m = 4, and use the approximate expression 
given in ( fTTI l. Using equation ( |25l l we get 



(a-l)ff//(e) w [To '^,exp(-£, 



--c,ep/£) + ^ exp(-£c,,r/e) + 2Ho\ 

(30) 

A comparison of numerical calculations of the effective 
time teff with equation (l30l l for values of a and m that satisfy 
a + OT = 4 is shown in figure |4] The effective time (multiplied 
by the speed of light) is shown in units of Mpc to show the 
effective distance from which these CRs can reach us. 
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Fig. 4. — Comparison of direct numerical calculations of the effective CR 
production time, with the analytic approximation 130) for a + m =4. The 
effective time is multiplied by a — 1 to show the similarity of the energy de- 
pendence and to compare with the analytic expression. The time is shown in 
units of Mpc. 

In fact, for energies e > lO'^eV, using equation ( |30] l, gives a 
good approximation (to better than 20%) to the effective time 
for values of m and a in the range 2 < a < 3 and < m < 3 
as shown in figure |5] 

3. ENERGY PRODUCTION RATE OF 
10'''eV < £ < lO^'eV UHECRS 

In this section we discuss the energy production rate of 
UHECRs required to explain the observed CR flux. Assum- 
ing that the CRs are mostly protons, we can use equation ( l30l l 
to estimate the effective generation time of CRs in the energy 
range e > lO'^eV. The energy production of CRs per loga- 
rithmic energy interval can be approximated by: 

e'fi(e) « 

(a-l) [r,^'^exp(-£c,fp/£) + %Jrexp(-£c,^/e) + 2//o] £^^, 

(31) 

where dn/de = J/(4ttc) is the measured CR particle density 
per unit energy. Assuming that the cosmic ray spectrum index 
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Fig. 5. — Comparison of direct numerical calc ulati ons of the effective CR 
production time, with the analytic approximation )30t for limited values of m 
and a. The effective time is multiplied by a - 1 to show the similarity of the 
energy dependence and to compare with the analytic expression. The time is 
shown in units of Mpc. 

is in the range 2 < a < 3 and that < m < 3 [see Eq. (|2|i], this 
equation is accurate to better than 20% (as explained below, 
systematic and statistic experimental errors lead to a larger 
uncertainty in Q). The energy generation implied by the latest 
data of Auger and Hires for e > lO'^^eV is shown in figure|6] 
As can be seen in figure |6] the energy generation rate in the 
range e > 10''*-^eV is e^Q(e) « lO"*^'^ - lO'*'^ erg Mpc"^ yr"'. 
It should be noted that some of the CRs in this energy range 
might be of Galactic origin (especially in the lower energy 
part of this range) and some of them may not be protons. In 
this respect, the energy production shown in the figure can 
be considered an upper limit to the energy production in XG 
UHECR protons. In addition, at energies e > lO^'^eV, the lo- 
cal CR source density, which may be differ ent from the aver- 
age one, ma y affect the measur ed flux (e.g.'Giler et al."1980'; 
iBahcall & Waxman 2000; Waxman 1995; Berezinsky 2008). 

The statistic and systematic errors in the generation rate es- 
timate are of the order of tens of percents. The main source 
for the systematic error is a 20% systematic uncertainty in 
the energy calibration of the CR particle energies, which im- 
plies a ^ 30% uncertainty in the measured flux. In addition, 
all experiments have statistic errors in the particles' energy es- 
timates. As shown in ap pendix § Ia1 the flux is overestimated 
by a factor of [Eq. (lA4l i1 

/stat « [1 +0.5(/3- l)(/3-2)aL] , (32) 

were cr^j^j = ((Asstat/e)^) and (3 = dlo^J/dlo ^s. For a value 
of (Tstat ~ 0.2 typical to these experiments (e.g. [Abraham et al] 
llOOS) this implies an eiTor of 30-40% at peak (e ~ lO^^eV). 
If higher precision is sought, an energy dependent estimate of 
(Tstat at these energies is essential. 

To illustrate the effect of these uncertainties, we show in fig- 
ure |7] the energy generation rate implied by the data of Auger 
and Hires with the Auger energy scale shifted by 20% (within 
the reported systematic error, and as required in order that the 
two experiments agree), and correcting for the statistic smear- 
ing using Eq.(l32ll with CTstat = 0.2 for both experiments. By a 
fractional shift, (Ae/£)sys, in the absolute energy scale of an 
experiment, we mean that particles with a measured energy e 
have in fact an energy of e + (Ae/e)syse. 

Another source of uncertainty is the unknown value of the 
spectral index a. The spectral index can in principle be di- 



rectly read off from figures |6] and [7] In practice, this is prob- 
lematic due to the unknown Galactic and heavier nuclei con- 
tribution at energies e < lO'^'^eV and due to the possible ef- 
fect of the unknown local CR source density on the flux at 
particle energies e > lO^^eV, as explained above. This leaves 
a limited range of energies where the observed spectrum can 
be safely used to infer the generation spectrum, 10'^ ''eV < 
e < lO^^eV, where the limited statistics and the unknown de- 
pendence of (Ae/£)sys and (Ae/e)stat on e does not allow an 
accurate determination of a. The spectral index cannot be 
much smaller than a « 2, as it would predict too many events 
at the highest energy and require a significant Galactic contri- 
bution at energies e > lO'^^. The spectral index is conserva- 
tively limited from above by a < 2.7, since larger values will 
produce a flux at lower energies, e ^ 10'^ eV, that exceeds the 
observed flux (see more detailed discussion in § 14.31 1. 

A simple analytic estimate of the XG energy generation 
rate can be made as follows. In the energy range lO'^^eV < 
e < lO'^^eV, slightly below the threshold for pion produc- 
tion by interaction with CMB photons, the energy loss of 
protons is dominated by pair production and the effective 
generation time is approximately given by t^jf « {tqIj, + 

2//())"' « 2.5 X 10^ yr. The energy density of CRs per log- 
arithmic particle energy at e = 10 '^'^eV is roughly e^dn/de « 
j^Q-20.6 gj-g (;jjj-3 xjius we can approximate the energy pro- 
duction at e < lO'^' ^'eV as: 

e^G(e < 10'''''eV) « {a- l)e^dn/de\io»-«^y x (t^i^^ + 2//o) 

« 0.3 X 10^^(a- 1) 33 erg Mpc yr . 

\ 10 ^"-^ erg cm / 

(33) 

Note that there have been claims in the literature that 
the production rate is consid erably higher or that the un- 
certainty is much larger (e.g. Farrar& Gruzinov 2008). The 
source of the confusion is that different authors estimate 
the production rate of XG CRs at different CR ener- 
gies, e.g. at £ > lO'^eV (as in this ms) vs. e ^ 
10^^ (e.g. iBerezinsiEvI 120081) . In particular, the genera- 
tion rate p er logarithmic p article energy in the model pre- 
sented by .Berezinskvl (1200 8.) at energies e lO'^^eV and 
e - lO'^eV is e^e(e)|io„.6eV « 1-3 x lO^* erg Mpc'^ yr and 
£^C(£)lioiseV ~ 17 X 10"^ erg Mpc"-^ yr respectively. The 
energy generation rate at e > lO'^eV is in agreement with 
the generation rate estimates presented here [given that in 
Berezinsky's model a = 2.7 and (AE/E)sys « 40% is required 
for the Auger energy calibration, in order that the Auger flux 
agrees with the model at e ^ 10''' eV, see § 14.311 . 

4. TRANSITION FROM GALACTIC TO EXTRA 
GALACTIC CRS 

In this section we examine the CR spectrum down to lower 
energies and discuss different models that may explain the 
smoothness of the spectrum from the knee to the highest ob- 
served energies. First we review in § 14. II the argument that 
a transition between two different CR sources should be ob- 
served as a flattening in the spectrum, which implies that the 
transition between Galactic and XG sources is likely to oc- 
cur at the 'ankle'. We demonstrate that a model in which the 
transition is at e ^ lO'^eV involves fine tuning. In § 14.21 we 
present two examples of specific models which are consistent 
with the data (figure|9]and fi gurefTOll. Finally w e show in § 14.31 
that the model presented in ( lBerezinskvl2008l) is somewhat in 
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Fig. 6. — The energy generation rate as measured by Auge r an d Hires 
assuming that the CRs are pu rely protons. We used equation BOt setting 
Q — 1 = 1. For < HI < 3, Eq. (30) is accurate to better than 20%. Statistical 
and systematic errors in the experimental determination of event energies lead 
to ~ 50% en'ors in the flux at the highest energies. For different values of a, 
the spectrum should be multiplied by an energy independent factor (a — 1). 
The absolute energy scales of the Auger and Hires data where not altered in 
this figure. 
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Auger (AE/E)^^^=0.2, 0^,^=0.2 
- HiRes (AE/E) =0.029 ,o = 



This requires fine tuning in both the spec tra and amplitudes of 
the two contributions (e.g. lHillaslll984l) . 

A class of su ch models pursued by some authors 
dBerezinskvl l2008l and references therein) suggests that the 
XG CRs are dominant down to much lower energies, e ~ 
lO'^eV. According to these models the flattening (ankle) is 
due to the flattening ('dip') in the CR energy loss time as a 
function of energy. In these models the transition from Galac- 
tic to XG sources has to occur without an observed flattening, 
which requires fine tuning between the parameters of the XG 
and Galactic components. To de monstrate thi s fine t uning we 
consider the model presented by iBerezinskvl (120081 figure 7 
left). A scanned version of this figure is reproduced in figure 
[8] On it, we added two slightly modified models (the blue and 
red thin full lines) in which the Galactic and XG components 
were enhanced by a factor of 3 respectively. As can be seen, 
there is a clear flattening that is expected in both cases at the 
transition region. Such flattening is not observed. 
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Fig. 7. — The energy generation rate as measured by Auger and Hires 
assuming that the CRs are purely protons. Here the absolute energy scales of 
the Auger and Hires data where altered as shown (in order that they agree at 
e ~ 10"eV, see 8 14.21 . In this figure we took int o ac count the statistic errors 
in the particle energy estimates according to eq. (32). 

disagreement with the latest published data from the Auger 
Experiment. 

4. 1 . The fine-tuning argument for a Galactic-XG transition 
at the 'ankle ' 

It is established that CRs at the lowest energies, e ^ 1 GeV, 
are of Galactic origin and it is likely that the highest energy, 
e > 10'^-^ particles are XG. A basic feature that is expected 
at any transition between two different sources of CRs is a 
flattening of the spectrum. This is expected since in order that 
the flux from the source of the higher energy particles exceed 
the flux from the source of particles at lower energies, it's 
spectrum must be flatter It is quite remarkable that there is 
only one observed flattening in the CR spectrum throughout 
the span of 10 orders of magnitudes of particle energies. This 
flattening, the so called 'ankle', is observed at an energy of 
roughly e ^ 10"*^eV. In any model where the transition is at 
lower energies, it must occur without an observed flattening. 



Fig. 8. — Demonstration of the fine tuning problem in models where the 
transition from Galactic to XG C Rs occurs without a flattening. The model 
presented in Figure 7 (left) from (Berezinsky 2008) is reproduced. The green 
asterisks, pluses and filled line are reproductions of the Galactic, XG and total 
fluxes respectively. The red (blue) filled line is a slight modification of the 
model where the Galactic (XG) contribution was enhanced by a factor of 3. 
As can be seen in both cases, a clear flattening that is not observed is expected 
in both cases. 

It is important to stress that the amplitudes and spectra of 
the two contributions in this model are affected by many phys- 
ically unrelated factors. In particular, these have to be differ- 
ent types of accelerators (the Galactic ones limited in energy 
and the XG capable of accelerating UHECRs). It is challeng- 
ing to offer a physical reason for the spectra to adjust in such 
a way as to erase the evidence for the transition. 

4.2. Two examples of models with a transition at the 'ankle' 

In figures |9]and[T0]we present two examples of phenomeno- 
logical models that include a Galactic and XG contribution 
that reproduce the observed spectra to a satisfying level. In 
figure |9] the Galactic contribution is a pure power law. The 
parameters are: 

7c = 1-7 X 10"^2(e/10"* ''eV)"^-2 m'^sec"' sr"' eV"', (34) 

e^Qxcie) = 0.45 x 10-^(e:/10'^-*eV)-''-2 erg Mpc"^ yr"' . 

(35) 
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In figure [To] the Galactic contribution is a power law with an 
exponential cutoff at Scf = 10'^ ''eV. The parameters are: 

sr"' eV" 
(36) 



7c= 1-7 X 10-^'-(e/10''*'eVr'-^exp(-e/£,./) m-^sec 



£ Qxcie) = 0.55 x 10-'*^(£/10'''"eVr"-^ erg Mpc"' yf 



(37) 
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Fig. 9. — A mo del of UHECRs with a power law Galactic contribution 
[Eqs. 1341 and )35H . The expected flu x is calculated numerically (full black 
line). The expected flux using Eq. )20t is shown for comparison for energies 
e > 10* ^e V (full red line). Th e measured flux ab ove the 'knee' a s measured 
by Akeno (Nagano & Watson"'200()), Fly's Eye (Bird et al."1994!, data taken 
from Nagano & Watson 2000) and HiRes (Abbasi et al. 2008) is presented. 
The absolute energy calibration of the various experiments shown were ad- 
justed in order that they agree at e ~ 10"eV. The fractional shifts in the 
absolute energy scale, AE/E, are within the published systematic en'ors. 
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Fig. 10. — Same as figurellOl with the Galactic spectrum apowerlaw 
an exponential cutoff, Eqs. )34t and j35t . 
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We make the following observations regarding these fits. 
These models agree with the data for energies e > lO'^eV 
to about 20%. There are many factors that affect the fluxes 
and are not taken into account: Contribution from heavier el- 
ements, diffusion in the intergalactic magnetic fields, non triv- 
ial contribution from nearby sources etc. We believe that with- 
out a better understanding of these systematics, using simple 
phenomenological models is the best one can do. On the one 



hand, it is encouraging that simple models can relatively well 
explain the data on such a wide energy range. On the other 
hand, it is clear that from studying the CR spectra it is very 
difficult to distinguish between the Galactic and XG contribu- 
tions. Thus, it is difficult to use these measurements to learn 
about the different unknown factors. The fact that two ba- 
sically different models agree with similar success is a good 
illustration of this problem. 

4.3. Can we already rule out some models based on the 
measured all particle spectrum? 

In order to use the all particle spectrum to rule out models, 
these models must have clear predictions for the spectrum. As 
the chemical composition of the CRs is unknown and since 
there can be a Galactic contribution to high energies (see § |4|i 
the all particle spectrum by itself gives us limited information 
on the detailed XG UHECR ge neration. 

The so called 'dip' models (iBerezinskyI l2008l and refer- 
ences therein), which assume a pure proton composition with 
no Galactic contribution at high energies, have a rather clear 
prediction for what the all particle spectrum should be. What 
makes these models interesting, is the claimed agreement 
between the observed spectrum and the predicted spectrum 
bas ed on this assump tion. Here we compare the model given 
in ( Berezin skyI l2008l) with the latest published data of the 
Auger experiment. This model has a CR energy production 
rate of £^Q{e) = 16.7 x 10'*4(e/10i**eV)-"-^ erg Mpc"^ yr"' 
(this corresponds to a total CR emissivity of Lq = 3.7 x 
10"*^ erg Mpc""* yr"' for particle energies e > 1 GeV, for a 



-2.7 



at 



spectrum Q oc e at 1 GeV < e < 10 eV and Qcc e 
e > 10"^eV). The spectral index and amplitude of this model 
were obtained by fitting the data from several experiments in 
the energy range e > 10'**eV. As claimed by the authors, 
the absolute energy calibration of the different experiments 
is fixed by the position of the 'dip' in the spectrum and thus 
the spectral index and the amplitude can not change by much. 

In figure [TT| the cumulative number of events above differ- 
ent energies measured by the Auger experiment (green aster- 
isks) is compared to the expected number according to this 
model (red line). As can be see n, the model over predicts the 
number of events (IStanevll2008h . To achieve agreement at en- 
ergies e ^ 10'''eV it is necessary to assume that the Auger ex- 
periment systematically underestimates the particle energies 
by about 40%. The spectrum that would result if the ener- 
gies of the particles were shifted by 40% is shown on the 
figure (blue dots). It seems that the shifted spectrum at the 
highest energies is flatter than expected by the model. For 
comparison, a model with a generation spectrum of Q (x e~^'* 
is plotted (red dashed line). It can be seen that even for this 
harder generation spectrum, which is too hard to be consistent 
with the observed spectrum at e lO'^eV, the resulting flux 
seems too soft to be consistent with the observed spectrum 
above e ~ 10'*^ "^eV. We stress that in order that this disagree- 
ment may be confidently used to rule this model out, the trend 
must be confirmed with higher statistics and the systematics 
of the experiment must be better understood. In addition, an 
enhancement of the flux at the highest energies, e > 10^'^eV, 
may be a consequence of an e nhancem ent in the local den- 
sity of CR sourc es (e. g.lGileret al. 1980: lBahcaU & WaxmanI 
I2000t IWaxman |[T99llBerezinskyii2008i) . 

5. SUMMARY AND CONCLUSIONS 

In this paper we derived simple analytic tools for the analy- 
sis of the UHECR flux and used them to estimate the UHECR 
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Fig. 1 1. — Cumulative number of events above different energies mea- 
sured by the Auger experiment (green asterisks) compared to the expected 
number according to iBerezinslcy 2008) (red line). The Auger spectrum with 
the particles shifted by 40% is shown on the figure (blue dots). 



generation rate. 

We first derived in section §|2] approximate analytic expres- 
sions for the effective generation time feff (e) = dn/d£/ Q(e,z = 

0) of CRs, assuming the flux is dominated by protons coming 
from XG sources. We used the fact that at high energies the 
energy loss time r is much shorter than the Hubble time //q ' 
to derive approximate expressions for the effective generation 
time, Eqs. (fT2] |. (V% . and (l20l l. by expanding it to first order 
in the small parameter Hqt. We then obtained a phenomeno- 
logical analytic approximation to the loss time of p rotons due 
to pai r production and pion production, Eq. ( |25] l, (IWaxman I 
Il995h . and used it to derive a simple analytic expression for 
the effective generation time, Eq. (|30] |. We showed that these 
expressions are in good agreement with detailed numerical 
calculations (figs. [T][3]|4]and|5]l. In particular, for 2 < a < 3 
and < OT < 3 [see Eq. dU], the simple expression given 
in Eq. ( [30l l agrees with numerical calculations to better than 
20% at £ > lO'^eV as shown in figure|5] 

Next, we estimated in §[3]the energy production rate of cos- 
mic rays with energies e > lO'^'^'eV using the latest reported 
measurements of the Auger experiment [figs |6] and |7] and Eq. 
^]. We showed that it is roughly energy and model inde- 
pendent and equal to £^2(£)|e>ioi95eV ~ 0.3-0.6 x 10'^'*(a- 

1) erg Mpc""* yr"'. This is consist ent with earlier result s de- 
rived by Waxmanl (|1995|) and Bahcall& WaxmanI (l2003h . and 
also by llerezinsky (200j), who obtains £^2(£) f« 0.75 x 

lO^(a-l)(e/lO^'^-^eV)-°'^^ erg Mpc"^ yr"' for a = 2.7. The 
slightly higher normalization obtained in the latter analysis 
is due mainly to the assumption, required by that analysis, 
that (Ae/£)sys — 40% for Auger (see § I4.3l l. a systematic er- 
ror which is larger than that quoted by the experiment (note, 
that the considerably higher CR energy production rate some- 
times quoted from Berezinskv 2008, refers to lower particle 
energies). 

The statistic and systematic errors in the generation rate es- 
timate are of the order of tens of percents. The main source 
for the systematic error is a ^ 20% systematic uncertainty 
in the energy calibration of the CR particle energies, which 
implies a ~ 30% uncertainty in the measured flux. Another 
source of uncertainty is the unknown value of the spectral in- 
dex a. There is a limited range of energies where the observed 



spectrum can be safely used to infer the generation spectrum, 
jQi9.5gY < £ < 10^°eV, where the limited statistics and the 
unknown dependence of (A£/£)sys and (A£/£)stat on e does 
not allow an accurate determination of a. a is roughly lim- 
ited to the range 2<a< 2.7 (see discussion in §|3]and § 14.3b . 

In section § |4] we showed that the measured flux of CRs 
at energies e > lO'^eV is consistent with simple models that 
include a simple Galactic contribution and a transition be- 
tween the Galactic and XG dominated regions at the 'an kle', 
£ ^ 10'''eV (figs |9] and [TOli. As explained in section § 14.11 
and illustrated in figure [8] models in which the transition be- 
tween Galactic and XG sources occurs at energies consider- 
ably below the ankle require a fine tuning in the amplitudes of 
the Galactic and XG contributions in order that the expected 
flattening in the transition will not be observed. Moreover, 
such models require a steep spectrum, a 2± 2.7, which is dis- 
favored by the data (see fig. fTTT l. In fact, even harder gener- 
ation spectrum, a = 2.4 , which is too hard to be consistent 
with the observed spectrum at £ ~ lO'^eV, seems too soft to 
be consistent with the observed spectrum above e ^ lO'^ ^'eV 
(fig. [TTT i. Nevertheless, it should be noted that in order that 
this disagreement may be confidently used to rule out models 
with Galactic to XG transition well below the ankle, the trend 
must be confirmed with higher statistics and the systematics 
of the experiment must be better understood. In addition, an 
enhancement of the flux at the highest energies, e > 10^"eV, 
may be a consequence of an enhancement in the local den- 
sity of CR sources (e. g.'Giler et al.l'1980t lBahcall & Waxmanl 
2000; WaxmanllI995l h Berezinsky 2008). 

To conclude, we summarize what we believe the main con- 
clusions that can be drawn from the analysis of the all particle 
spectrum are. 

1 . The detection of the GZK cutoff (Bahcallcfe Waxmanl 
|2003t lAbMsLet al,| |2008l [Abraham et al. 2008) sup- 
ports the claim that UHECRs with energies e ^ 10'''eV 
are mainly protons arriving from XG sources. We note 
that the fact that we obtained a reasonable, smooth, CR 
generation spectrum (figs |6] and |7]i after correcting for 
the particle energy losses (which are strongly energy 
dependent) is consistent with the presence of the GZK 
cutoff. 

2. Assuming UHECRs are XG protons, the energy gener- 
ation rate is e^Q(e > lO'^-^eV) ~ 0.3-0.6 x lO^^V" 
1) erg Mpc"' yr"' [see figs |6] and Hand Eq. ^]. 

3. The fact that the only flattening observed in the spec- 
trum is at the 'ankle' strongly suggests that the tran- 
sition from Galactic to XG CRs occurs at that energy 
scale (see discussion in § 14. Il l, lO'^eV. 

4. The data in the energy range e > lO'^eV is consistent 
with simple models including Galactic and XG contri- 
butions (the latter with a flat energy generation spectral 
index, a 2± 2) where the transition occurs at the 'ankle' 

(figsisnoj. 

5. The all particle spectrum alone cannot be used to dis- 
cern the detailed spectral shapes of the Galactic and XG 
contributions (see § I4.2| |. There is a limited range of en- 
ergies where the observed spectrum can be safely used 
to infer the XG generation spectrum, lO'^^eV < £ ^ 
10^"eV, where the limited statistics and the unknown 
dependence of (A£/£)sys and (A£/£)stat on e allow a 
broad range for the spectral index, 2 < a < 2.7. 
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APPENDIX 

THE EFFECT OF STATISTIC ERRORS IN THE MEASUREMENTS OF INDIVIDUAL PARTICLE ENERGIES 

ON THE DERIVED FLUX. 

Suppose that each particle with energy eo has a probability P(^-^)^ to be measured to carry an energy e in the interval 



e,e + de. The number of particles measured in the interval [ei , £2] is 



N{ei <e <e2) = 



e, ^0 



deo—{eo)P( ) — , 

de eo eo 



(Al) 



where dN /de is the number of particles per unit energy that reached the detector during the measurement. By changing variables 
[e, eo] to [e, jc = we can rewrite this expression as (after some straight forward algebra) 



de j dx — ( ) P(x)dx. 

_i de l+x l+x 



(A2) 



Assuming that the distribution of arriving particles can be described locally (in energy) as a power law distribution dN /de = 
Ae~^, where the value of (3 can be derived from the derivatives of the distribution in the energy regime under consideration, 
f3 = dlogiN)/dlog{e) we find [using Eq. dATH 



N(ei <e<e2)=A / dee~^ / dx(l +xf-^P(x)dx = Noiei < e < £2) J dx(l +xy-'P(x)dx, 



(A3) 



where A^o(£i < £ < £2) is the true number of particles that reached the detectors in the range £; < £ < £2. Assuming asm = (x~'^ <C 
1 and that the energies are calibrated so that (£) = £0, we get 



N(e, < e < £2) = A?o(ei < e < £2)(1 + <^L + Oiat,J). 



(A4) 



The value of (3 can be written in terms of the production and propagation as /3 = a + /3prop , where a = -d log(2) /d log e and 

(3prop = -dlogiteff)/dloge. (A5) 

For propagation of protons, (3prop reaches approximately 3 at £ ^ 10^°. For 2 < a < 3, /? can be as large as 5 to 6. For a value 
of (Tstat ~ 0.2 typical to these experiments (e.g. [Abraham et aT]|2008l) this implies an error of 30-40% at peak (e ^ 10^°eV). If 
higher precision in the flux estimate is required, an estimate of a at these energies is essential. 
At lower energies e < lO'^^, we measure /3 w 3 directly and the errors are of the order of a few percents. 
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